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1. Introduction 

In light of the apparent anisotropy of the cosmic background radiation map (CMB) 
obtained by space missions, the investigation of anisotropic cosmological models seems 
to gradually regain scientific interest almost forty years after its peak in seventies 
PP,[2]- These models are described by a three dimensional homogeneity group and a 
trivial or reduced isotropy group as compared with Friedman-Robertson- Walker (FRW) 
models. Of particular interest are different scenarios of the early epoch of the universe 
including also quantum effects. One theoretical paradigm that is adequate for describing 
such situations is quantum field theory in curved spacetimes (QFT in CST). This 
theory pertains mainly to mathematical physics and aims at a rigorous description 
of physical phenomena starting from an axiomatic basis. The way from the axioms 
to explicit results readily applicable at the observational level usually requires a rich 
supply of mathematical methods. In the context of (anisotropic) cosmological models, 
harmonic analysis, and related, Fourier analysis with respect to the isometry group of 
the spacetime under consideration are powerful tools for obtaining explicit results. 

There is a prevalent tendency among different scientific disciplines to specialize to 
their own aims and perspectives and thus diverge from each other. This results in an 
isolation of scientific results inside a limited community without a wide access from 
outside. Perhaps one of such fields is the mathematical theory of harmonic analysis. 
The Levi decomposition effectively breaks apart the general harmonic analysis to those 
for solvable and semisimple groups separately. Both branches have been investigated 
in great generality. The Kirillov orbit theory for solvable groups, along with such 
developments as the Currey theory for exponential solvable Lie groups, has given 
methods for obtaining explicit results for arbitrary dimensions. On the other hand the 
theory developed by Harish-Chandra, Helgason and others provides a very deep insight 
into general semisimple homogeneous spaces and related structures. At the same time, 
abstract harmonic analysis for locally compact groups, and the Mackey machine based 
upon it, yielded remarkably in describing general principles and phenomena. However, 
when a theoretical cosmologist wants to perform a mode decomposition of some physical 
field on a homogeneous spacetime it is almost of no use to him to know that, for instance, 
a Borel measure exists, or that it can be computed up to equivalence for N dimensions. 
What a cosmologist really needs is an explicit description with formulas which can be 
used without expert knowledge in harmonic analysis, as it is available for the traditional 
Abelian group M n . But as far as we were able to see, results of that kind do not seem 
to be available in the literature; possibly because they are usually not in the focus of 
interest of mathematicians. It may be regarded as one of the duties of mathematical 
physics to provide bridges between the increasingly diverging interests of mathematics 
and physics. With this in mind we take to the task of giving an explicit description of 
harmonic analysis of a number of groups of cosmological interest. 

As stated above, in quantum field theory on cosmological spacetimes one is often 
interested in Fourier analysis with respect to the symmetries of the underlying geometry. 
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In the framework of homogeneous cosmologies, M = IxE where I is an open interval, 
and £ is a smooth three dimensional manifold with a Riemannian metric h on it. The 
spacetime metric g on M takes the line element form 

ds 2 = dt 2 — aij^ds^ds^ 

with a positive definite matrix aij(t) depending smoothly on the "time coordinate" 
tel. (£, h) is a homogeneous space with respect to a Lie group of isometries G. 
The variety of all such homogeneous spaces which arise in cosmology can be found in 
several works including [3], [I], [5]. We have analyzed the geometrical setup of general 
linear hyperbolic fields on cosmological spacetimes in [6]. In particular we have seen 
that in nearly all cases one deals with a semidirect product group G = X x O where 
O is either of SO (3), SO (2), {1}, and the role of £ is played by Bianchi I-IX groups 
Bi(N) (N = I , II , IX) and their quotients Bi(N)/T by discrete normal subgroups 
T. Such structures are called semidirect homogeneous spaces, and a few important 
results have been obtained in [6] in this generality using abstract harmonic analysis. 
We moreover have seen how tightly harmonic and spectral analysis is related to mode 
decomposition. The next step towards physics will be to describe harmonic analysis 
of all those possible semidirect homogeneous spaces explicitly. The spaces of maximal 
symmetry with O = SO(3) are the FRW spaces, which are described by isometry groups 
50(4), E(3) or £0 + (l,3). The spaces with one rotational symmetry are described by 
O = S0(2) and are called LRS (locally rotationally symmetric) spaces. And finally the 
purely homogeneous spaces are given by trivial isotropy groups O = {1}. The isometry 
groups of FRW spaces are classical groups and their harmonic analysis is also a classical 
subject. For purely homogeneous spaces (otherwise called Bianchi spaces) this is known 
partially. The Bianchi I group is the additive group M 3 of which harmonic analysis is 
textbook standard. The Bianchi II group is the famous Heisenberg group of dimension 
2 + 1, which is well studied, and its harmonic analysis can be found in [7], [8], [9]. The 
Bianchi III group is the ax+b group in 2+1 dimensions whose harmonic analysis is known 
as well [9]. The Bianchi VIII group is the universal covering group SL(2, M), which is a 
notorious non-linear group. Its harmonic analysis can be found in pi)]. The last group, 
Bianchi IX, is simply SU(2) which is again classical. Little is known about the Bianchi 
IV- VII groups beyond the structure of their Lie algebras which are semidirect products 
of Abelian algebras M 2 and K.. In fact, although there are principally no obstacles on the 
way of their investigation, we were not able to locate any explicit description of their 
harmonic analysis in the literature; we also asked some prominent experts in the field, 
and none was able to point to such references. Even less is known about the semidirect 
products of Bianchi groups with 50(2) describing LRS models. We stress again that 
there is no obstacle to applying the Mackey machine and perform all calculations, but 
it seems that this has not been done so far. Because the discrete subgroups T can be 
readily found from the group structure [H] , once having control over harmonic analysis 
of the group G it is not hard to reduce it to the quotient G/T, but it again needs to 
be done somewhere. We have chosen to start with harmonic analysis of the solvable 
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Bianchi I- VII groups in a uniform manner deferring the remaining structures to the 
future. One could argue that this might be done straightforwardly by the exponential 
solvable methods of [12], but if one actually starts to do that (what we indeed did) one 
needs to construct an enormous amount of spaces, dual spaces and intersections, which 
are designed to handle arbitrary groups, and seem to be too bulky to be performed by 
hand. Therefore we preferred the original Mackey construction. 

The harmonic analytical Fourier transform is only half the way to applications. In 
fact what is really useful in concrete computations is the spectral Fourier transform 
given by the eigenfunction expansion of a left invariant Laplace operator A acting 
on the sections in the subbundle of T over S. In [6 J we have seen how the mode 
decomposition of linear fields can be performed knowing the eigenfunction expansion of 
A. We moreover have seen how in principle one can identify the harmonic analytical 
Fourier transform with the spectral Fourier transform of A to translate the general 
results on semidirect spaces into the language of concrete calculations. For this purpose 
one needs to find the spectral theory of A explicitly for all possible £ and invariant 
Riemannian metrics on them. This is again something that can hardly be found in 
the literature, although spectral analysis in Riemannian spaces is huge and very well 
developed a subject in mathematics. In particular, one needs to know the spectrum 
and a complete system of eigenfunctions of A explicitly. In general, the eigenfunction 
problem of A is a vector valued elliptic partial differential equation on a manifold 
without boundary, which is difficult to compute even numerically. If this equation 
admits separation of variables so that the eigenfunctions are given by combinations of 
functions of one variable subject to ordinary differential equations then we can consider 
these eigenfunctions as given explicitly in terms of special functions. In this work we 
will give such an explicit description of the spectrum and eigenfunctions of A in terms 
of an arbitrary left invariant Riemannian metric on £ for the line bundle over Bianchi 
I- VII groups. For arbitrary bundle dimension this is much more complicated. There 
is a bit of hope to obtain explicit solutions by transforming the original vector valued 
eigenfunction equation on the manifold to a scalar elliptic eigenfunction equation with 
constraints on the holonomy bundle of the linear connection associated with the given 
fiber metric. This is a non-trivial task to which we hope to return in the future. 

We summarize the content of the current exposition as follows. First the Bianchi 
I- VII groups are realized as semidirect products of R 2 and K. and the main group 
properties are explicitly computed, such as the multiplication laws, exponential maps, 
modular functions and adjoint representations. Then the dual spaces of the groups are 
constructed, i.e., the equivalence classes of unitary irreducible representations. This is 
done by means of the Mackey machine. Next a look is given at the co-adjoint orbits 
of the groups in the sense of the Kirillov theory, and it is described explicitly how the 
cross sections can be chosen to parameterize the dual space. Afterwards an explicit 
Plancherel formula is given for all these groups. Thereafter we turn to spectral analysis. 
The spectra and the eigenfunctions of A are found explicitly in terms of the chosen 
arbitrary left invariant Riemannian metric. Then it is shown that these eigenfunctions 



Harmonic analysis in Bianchi I- VII 



5 



are complete in L 2 and give rise to a conventional Fourier transform in sense of [5J. In 
the final part the applications in quantum field theory are discussed. First it is shown 
that Bianchi spacetimes are ideally adapted for the mode decomposition as given in 
[6J . Then this mode decomposition is demonstrated in the example of the Klein-Gordon 
field on Bianchi I- VII spacetimes. A number of interesting consequences are indicated 
including those for the quantum Klein-Gordon field, where results from [13] are used as 
well. 

2. Semidirect structure of Bianchi I-VII groups 

As a first step we will try to explicitly realize the solvable Bianchi II- VII groups (I is 
Abelian and will serve as a starting point in the analysis of others) as semidirect products 
of Abelian subgroups. A classification of solvable real Lie algebras with respect to such 
products can be inferred from |14j . 

Semidirect products of Lie algebras and Lie groups. We start by recalling 
some definitions. Let a and b be Lie algebras, and let D(a) be the Lie algebra of 
derivations on a. Let further / : b i— y D(a) be a Lie algebra homomorphism. The 
semidirect product Lie algebra a x f b is the algebra modelled on a © b with the Lie 
bracket 

[(a, 6), (a', b')} = ([a, a'] + f(b)a' - /(&>, [b, 6']), (a, 6), (a', b')ea® b. 

Let, on the other hand, A and B be Lie groups, and F : B h-> Aut(A) a Lie group 
homomorphism (Aut(A) embedded into GL(A)). The semidirect product A x p B of 
groups A and B is defined as the Lie group modelled on the product manifold A x B 
with the multiplication 

(a, b)(a', b') = {aF{b)a', 66'), {a, b), (a', b') e A x B. 

Following the notation of [15], denote by F° : B i— >• Aut(a) the map B 3 b ^ d[F(b)\ e 
Aitt(a), where a is the Lie algebra of A. Then the derivative of this map, / = will 
be a Lie algebra homomorphism / : b i— > D(a) (b the Lie algebra of B), and the Lie 
algebra of the direct product Lie group A x F B is the direct product Lie algebra a x f b 

Bianchi I-VII groups as semidirect products. With this in mind, let us start 
with realizing Bianchi I-VII algebras as semidirect product algebras g = R 2 x f M with 
some Lie algebra homomorphism / : R i— > D(M. 2 ). This correspondence between Bianchi 
algebras and homomorphisms / can be obtained by combination of [16] and [H]. (Those 
uncomfortable with Russian may simply perform the semidirect product construction 
and check the commutation relations.) Namely, in each case f(r) — r ■ M, r 6 R, in a 
suitable basis, where M is a 2 x 2 matrix. The matrix M for each algebra is given in 
Table [1] below. 

The corresponding integral homomorphisms F° will be the exponentials F°(r) = 
e rM (note that the exponential map on the group R is given by the identity map). If a 
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(::) 


(o - q ) 

-1 < q < 1 


(i 1) 

p>0 



Table 1. The matrices M for Bianchi I- VII groups 
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III 


IV 




VI 


VIII 
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(r l) 


(•; :) 


U o\ 

\re r e r 1 


(: ."■) 


|V \ 

I e~ qr j 
-1 < q < 1 


/ e pr cos(r) — e pr sin(r) \ 
I e pr sin(r) e pr cos(r) J 
p>0 



Table 2. The matrices F(r) for Bianchi I- VII groups 



diffeomorphism is given locally by a linear coordinate map, x\ = A\xj with the matrix 
A, then its differential will be given by the same matrix A. Now that F°(r) = d[F(r)] 
and that F(r) are linear automorphisms, it follows that F(r) = e rM . Thus all Bianchi 
groups I- VII are given by semidirect products G = M 2 Xp K, where for each class the 
group homomorphism F:1h^ Aut(M. 2 ) is given as in Table |2] above. 

We appoint to use capital symbols X, Y, Z for Lie algebra coordinates and small 
symbols x, y, z for Lie group coordinates, but these may interfere in some calculations 
involving exponential maps. It follows that the group multiplication is 

{x,y,z){x',y',z') = {{x, y) + F{z){x', y'),z + z'), (x,y,z),(x',y',z') 6G = K 2 x f l. 

The exponential map. Finally we note that all 7 groups are exponential, and 
the exponential map is given as follows. Let (X, Y, Z) G g = M 2 x f R with (X, Y) E R 2 
and Zef. We use the Zassenhaus formula 

exp(yl + B) = exp(A) exp(B) exp(C 2 ) exp(C 3 )..., 

where the coefficients C m are homogeneous Lie algebra elements composed of nested 
commutators of order m. We will use the convenient method of obtaining C m recursively 
as given in [T7]. If we set A = (X, Y, 0) and B = (0, 0, Z), we obtain 

[A B] = -f(Z)A. 

Now equating the homogeneous summands of any order of (4.7) and (4.8) of [17] , we 
obtain recursion formulas for C m which are bulky in general. However, trying an ansatz 
C m = a m (—f ) m ~ 1 (Z)A, a m G R, and checking directly for m = 2, one can easily prove 
it inductively, and find 

1 — m 



It remains to calculate 

exp(C 2 )exp(C 3 )... = exp (f^ ^-f) m ~^ Z ) A 

\m=l 
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If f(Z) is invertible for all Z then we write 

ml ml [m — \)l 

and obtain 

D(Z) = J2 ^(-f) m -\Z) = (-f)-\Z) (e-W - 1) - e-™ = 
A — ' ml 

m=l 

= f-\Z)(l-F(-Z))-F(-Z). (1) 

It is only for Bianchi II and III that f(Z) is degenerate, and for these two we can 
compute directly 

D{Z) = ^{-f) m -\Z) = -f(Z) for Bianchi II 

L — ' ml 2 

m=l 



oo 

1 — m 



and 

D(Z) = V ^^(-f) m - l (Z) = (1 - 2e- 1 )f(Z) for Bianchi III. 
^— ' m! 

m=l 

Thus we arrive at 

exp((X, Y, 0) + (0, 0, Z)) = exp((X, Y, 0)) exp((0, 0, Z)) exp(D(Z)(X, Y), 0). 

The exponential maps of M 2 and R are the identity maps, therefore 

(x,y,z)=exp((X,Y,Z)) = (X, Y, Z)(D(Z)(X, Y), 0) = ([1 + F(Z)D(Z)](X, F), Z), 

where F(Z) should be understood as F(exp(Z)). The matrices D(Z) appear somewhat 
bulky so we refrain from presenting them in a table. 

The adjoint representations Ad and ad. Let (g x , g y , g z ), (x,y, z) G G. Their 
conjugation (x',y',z') = (g X} g y ,g z )(x } y } z)(g X} g y} g z )- 1 is given by 

Or', y', z') = ((1 - F(z))(g x , g y ) + F(g z )(x, y),Z). 

The adjoint representation Ad is the differential of this map at the identity element 
(x, y, z) = (0, 0, 0), and so it is given by the matrix field Ad g , 



Ad g 



(F{g z ) -F'(0)(g x ,g y y 



\ 



1 

The adjoint representation of the Lie algebra is given by the matrbj|| 

ff(Z) -f'(0)(X,Y)A 

ad(x,Y,z) — 

y J 

% We use the general relation adx Y=[X,Y] for elements X, Y in a general Lie algebra. 
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The Haar measure and the modular function. The Haar measure on the Lie 
group is given by 

d(exp(X,Y,Z)) = j(X,Y,Z)dXdYdZ, 

where 

1 _ e - ad (X,Y,Z) 

j(X, Y,Z) = i) det , (X, Y, Z) G 0, 

ad(x,Y,z) 

and < [} 6 1 is an arbitrary constant. In group coordinates one can check that it is 
given by 

dg = t) det F(-g z )dg z dg y dg z , (g x ,g y ,g z ) G G. 

The groups are all non-compact, so there is no preferred normalization for the constant f). 
Later it will be determined as related to the chosen left invariant Riemannian metric on 
G. The modular function A(g) = det Ad~ 1 can be readily seen to be A(g) = det F(—g z ). 

This temporarily completes our task of analyzing the Bianchi I- VII groups as 
semidirect products. In the next section we will concentrate on their dual spaces. 

3. The irreducible representations of Bianchi I- VII groups 

In this section we will try to find the dual spaces of Bianchi I- VII groups using the 
Mackey procedure. Let us start with Bianchi I, which is simply the additive group M. 3 . 
Its dual group M 3 is homeomorphic to itself, M 3 = M 3 , and the irreducible 1-dimensional 
representations are given by 

= e i{ ^ } , x G M 3 , fcel 3 = M 3 , 

where we appoint to denote by {a, b} the usual Euclidean product of three-vectors 
a, b G M. 3 . These scalar functions £g can be viewed as unitary operator valued functions 
acting on the one complex dimensional Hilbert space C. 

The Mackey procedure for normal Abelian subgroups. We cite here the 
setup of the Mackey theory for groups with a normal Abelian subgroup as given in [9]. 
Let G be a locally compact group and N an Abelian normal subgroup. Then G acts on 
N by conjugation, and this induces an action of G on the dual group N defined by 

gv{n) = v(g~ 1 ng), g G G, v G N, n G N. 

For each v G N, we denote by G v the stabilizer of u, 

G v = {g G G: gv = z/}, 

which is a closed subgroup of G, and we denote by O v the orbit of v. 

a = {gv- g e G}. 

The action of G on N is said to be regular if some conditions are satisfied. To avoid 
presenting excessive information we only mention that if G is second countable (which is 
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true for a Lie group), then the condition for a regular action is equivalent to the following: 
for each v G N, the natural map gG u t— > gv from GjG v to O v is a homeomorphism. In 
our case N is a smooth manifold, and the group actions are all smooth, hence this map 
is not only a homeomorphism but even a diffeomorphism. The constructions become 
simpler under the assumption that G is a semidirect product of iV and the factor group 
H = G/N. We define the little group H v of v G N to be H v = G v fl H. Now we cite 
a beautiful theorem which appears as Theorem 6.42 in [9j and expresses the essence 
of the Mackey procedure. The functor Ind and the inducing construction are briefly 
introduced in the Appendix A. 

Theorem 1 (Folland,|9j) Suppose G = N tx H , where N is Abelian and G acts 
regularly on N. If v & N and p is an irreducible representation of H u , then Ind^ (i/p) is 
an irreducible representation ofG, and every irreducible representation ofG is equivalent 
to one of this form. Moreover, Ind G (yp) and Ind G (z/p') are equivalent if and only if 
v and v' belong to the same orbit, say v' = gv, and h i— > p{h) and h i-> p'(p _1 /ip) are 
equivalent representations of H v . 

Application to the Bianchi groups. It is easy to see that Bianchi groups II- VII 
satisfy the assumptions of the theorem. In this case N = M 2 and H = R, the dual of N 
is N = M? and is given by 

TV = {e i{ ^ } : x/keM 2 }, 

where we overload the notation by brackets {a, b} to denote the two dimensional 
Euclidean product of a, b G M 2 . Let «)y : R 2 4 G be the natural inclusion. The 
action of G on N is given by 

gfy{x) = ^ k (i N \g-h N (x)g)). 

All Bianchi solvable groups are homeomorphic to M 3 , and we may choose a global chart 
on them. In particular we choose one adapted to the semidirect structure M 2 Xpl 
presented in the previous section. Then the multiplication law in G is given by 

(x,y,z)(x',y', z') = ((x, y) + F(z)(x', y'),z + z ). 

The unit e G G is given by e = (0, 0, 0), and the inverse map by 

(x^^y 1 = (-F~ 1 (z)(x,y),-z). 

In particular, if (x, 0) = (x,y,0) G «at(M 2 ) and {g x ,9y,9z) G G, then 

(9x, 9 y , 9 z y 1 (x,y,0)(g x , g y , g z ) = {F' l {g z ){x, y), 0), 

that is, the conjugation map n >-)■ g~ l ng is given by (x,y) i-> y). Thus the 

action of G on iV is 
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where F- 1 ^) is the inverse transpose of the matrix F{g z ). This means that this action 
can be described by 

gk = F L (g z )k, g eG,keR 2 . 

Denote by V° C R 2 the eigenspace of M T corresponding to the eigenvalue (the null 
space). Then it will be also the joint eigenspace of the matrices F L (g z ) = e~ 9zMT 
corresponding to the eigenvalue 1 simultaneously for all g z G R. Let us write the 
stabilizer condition, 

e -^ MT k = k. 

Then the stabilizer G% and the little group will be 

Gk = MR 2 ) • Hi 

and 



R if k G V°, 
{0} else. 

Define the following space of irreducible representations of G: 

j=(V°x R)U(R 2 \V°). 
For each /j, G J the corresponding irreducible representation is given by 

TM = e'&ftj*"*, pt = (k, h) = (h, k 2 , k 3 ) 

if n e V° x R, and 

T, = T k = IndUe tfK} )^ = k 

if fi G R 2 \ V°. The orbit Of, is {k} if k G V° and F ± (R)k otherwise. As mentioned in 
the theorem, two representations /i, ji' G J are equivalent if and only if k and k' are on 
the same orbit, k = F ± (z)k', and the corresponding representations of H%. and Hy are 
equivalent when intertwined with the action of z. The first condition can be satisfied 
non-trivially if k,k' G R 2 \ V°, but then H~ k = Hy = {0}, and thus there exists only 
the trivial representation p — 1. Thus representations //,// G R 2 \ V° are equivalent if 
and only if they are on the same orbit. On the other hand, let //, // G V° x R such that 
= k', and the first condition is satisfied trivially. Then G^ = G, and G/G%. = {1}, 
so the action of 1 cannot intertwine inequivalent representations of Hy Thus // ~ // 
means = //. Therefore the dual space G of G will be 



_L 



G = (V° xl)U (R 2 \V°)/F 

The null spaces V°. Finally let us find the eigenspaces V° for different Bianchi 
groups. By a calculation of eigenvectors and eigenvalues of M we obtain 

Vj° = R 2 , V# = R © {0}, Vfjj = {0} © R, 
V? v = {0}, = {0}, Vyjj = {0}, 
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and 




{0}©R if g = 0, 
{0} else. 



Note that as always with solvable groups, the irreducible representations are either 
1-dimensional or infinite dimensional. 

To obtain explicit descriptions of the dual groups G for each Bianchi class we have 
to calculate the orbits Of. = F ± (M)k explicitly, which is done in the next section. 
Note that the entire construction could have been performed through the machinery of 
exponential solvable Lie groups developed in [TS] and [12], where the problem is treated 
exhaustively. In particular, it was shown that (as adapted to our terminology) there 
exists a cross section K , an algebraic submanifold of R 2 which crosses each generic orbit 
(i.e., an orbit of maximal dimension) exactly once, and thus parameterizes the infinite 
dimensional representations. Having explicitly calculated K we find G = (V° X R) U K. 
But the methods of [18] are extremely general and involve simple but lengthy algebraic 
calculations; this is why we have preferred the original topological Mackey constructions. 



4. Co-adjoint orbits of Bianchi II- VII groups 



The term co-adjoint orbits would probably suit better to the solvable Lie theoretical 
method of orbits as established by Kirillov and accomplished by Currey. At this point 
we deviate to present a little digression demonstrating the equivalence of that approach 
with that we have adopted. 

The Kirillov approach. The Lie algebra g = R 2 x / R of G is modelled on the 
vector space R 3 , and as such its dual space g' is again isomorphic to R 3 . We will fix 
this isomorphism by choosing the basis in g' dual to our adapted basis of g. With 



this identification the co-adjoint action of G on g' 
Ad* = Ad^r, 



is given by the matrix field 



Ad* 



For any I = (X*,Y*,Z* 



\(g x ,g y )F ± (g z )M T 
G g' its orbit 0\ is given by 





V 



O x = (F ± (R)(X*,F* 



)M T {X*,Y*) + Z*), 



and the space of orbits {0\} with the quotient topology induced from g' is homeomorphic 
to G with the Fell topology [9]. One can easily see that the orbits are of two types: 
those of (X*,Y*,Z*) with (X*,Y*) e V° or (X*,Y*) V°. The former are the so 
called degenerate orbits with dimension (singletons), and the latter are the generic 
orbits with maximal dimension 3. This is exactly the same result we obtained above by 
Mackey machine. 

The generic orbits and the cross sections. We will denote the range of a 
parameterized quantity Q(p) of a parameter p G P by Q(P). For instance, F J -(R) will 
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denote the range of the quantity F ± (r) when r runs over R. Here we will try to find 
the generic orbits F- L (R)A;o E G mentioned in the previous section and corresponding 
cross-sections K E R 2 . The latter will be algebraic manifolds composed of one or more 
connected components. In all cases V° is a subset of Lebesgue measure in R 2 . By 
the definition of the cross section K, the submanifold R 2 \ V° can be parameterized 
by a global chart k = k(k,r), (k,r) E Si x R, such that k(k,r) = F- L (r)k (k) and 
k (k) = k(k,0) E K. Under this diffeomorphism the Lebesgue measure dk becomes 
p(k,r)dkdr, where p(k,r) = \detd(k)/d(k,r)\. 

Now let us proceed to the determination of the orbits and the cross sections case 
by case. Figure 1 in Appendix B illustrates them qualitatively. 

II. We have 

F {r){k x , k y ) = {k x — rk y , k y ), 

hence the orbit through k E R 2 \ V° is F ± (K)(k x , k y ) = (R, k y ). The cross section can 
be chosen to be K — k (M) : ^ = R \ {0}, k (k) = (0, k). Indeed, any orbit (R, k y ) meets 
K exactly once at ko(k y ). Then 



p(k,r) = 
III. In this case 



. x dUk) dF ± {r) r .,. 



\k\ 



F' L {r)(k x , k y ) — (e r k x ,k y ), 

and the orbit through k G R 2 \ V° is F ± (R)(k x , k y ) = (sgn(k x )R + , k y ). Let Si = 
R x {—1,1}, k = (ki,k 2 ). The cross section is the image (— 1,R) U (1,R) of the map 
k (k) = (k 2 ,h). We find 

p(k, r) = e~ r . 

IV. For this group 

F {r)(k x , ky) — (e k x re k y , e k y ), 

and the orbits are complicated. We set Si = R +0 x { — 1,1}, k = (ki,k 2 ) and 
ko(k) = (k2,k 2 ki). That this is a cross section can be checked immediately. The 
measure density p is 

p(k,r)=e- 2r (l + k 1 ). 

V. Now 

F (r){k x , ky) — e (k x ,k y ), 

and the orbits are simply the incoming radial rays. Set Si = R/27rZ and k (k) = 
(cos(fc),sin(A;)). It follows that 

p(k, r) = e~ 2r . 

VI. For this group we consider only the case q ^ as q = is just the group III. 

F~ L {r)(k x , ky) = (e r k x , e qr k y ), 
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and the orbits are incoming polynomial curves if q < and hyperbolic curves if q > 0. 
For q < set & = R/2vrZ and k (k) = (cos(fc), sin(fc)). Then 

p(fc,r) = e-^-^ r (coB 2 (ife) - gsin 2 (A;)). 

For q > set £ = R +0 x {0, 1, 2, 3}, k = (k u k 2 ) and 

ko(k) ■ 



k 2 

— 1 \ / 1 



1 / Ui 



Thus 



p(k,r) = q k2 m°d2 e -(l- g )r^ 

VII. The co-adjoint action in this group is given by 

F ± (r)(k x , k y ) = e~ pr (k x cosr — k y sinr, k x sinr + k y cosr), 

and the orbits are incoming or outgoing spirals depending on whether p < or p > 0. 
We take ^ = (— e' rp ; —1] U [l,e' rp ) and fc Q (fc) = (fc, 0). Each orbit clearly intersects K 
exactly once. Finally 

p(k,r) = e~ 2pr \k\. 

Note that in all cases we have chosen K such that it possesses an involution 
ko(—k) = —ko(k), which will be useful in later constructions. Of course, these choices 
of cross sections are not unique, neither need they correspond to those suggested by 
Currey theory. In fact, one may make any other choice for convenience and calculate 
the corresponding measure density p precisely as we did. 



5. The explicit Plancherel formula for Bianchi II- VII groups 

We will obtain the Plancherel measure by extending the idea suggested in [9] for 
Heisenberg groups to all solvable Bianchi groups. Namely, we will exploit the Euclidean 
Parseval equality on the homeomorphic space M 3 . 

Introductory material. Before going to the solvable groups II- VII let us recall 
the well-known form of the Plancherel formula for the Abelian group R 3 . The Fourier 
transform of a function / G C^°(1R 3 ) is defined by 

/(*) = f dxe-^f(x), 

and the Plancherel formula is 

I dx\f(x)\ 2 = (2n) 3 [ alk\f{k)\ 2 . 

The Plancherel measure is simply dv(k) = (2ir) 3 dk, proportional to the Lebesgue 
measure on R 3 . 
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We start by noting that, being an algebraic (matrix) group, G is necessarily type 
I (Theorem 7.8 or 7.10 [9]), and the normal subgroup iV is unimodular and therefore 
in the kernel of the modular function A. It follows from Theorem 7.6 in [9] that 
the Mackey Borel structure on G is standard, and thereby due to Lemma 7.39 in 
[9] we have a measurable field of representations n p on p G G, such that tc p G p (or 
equivalently, we have a measurable choice of representatives of each equivalence class 
[n] G G). Henceforth we will speak of a representation 7r G G meaning the value of this 
field at a given point [ir] G G. As can be inferred from [T2j in the language of solvable Lie 
groups, only those irreducible representations corresponding to the generic orbits (i.e., 
orbits of maximal dimension) admit a non-zero Plancherel measure. Therefore only T M 
with n G K 2 \ V° ( generic representations) will play a role in the Fourier transform. We 
proceed to their construction as T k = Ind^e*^'^) following §6.1 in [9]. 

The Fourier transform at generic representations. For each k G R 2 \ V° 
the representation Hilbert space H k of ut = e^''' is H% = C. The homogeneous space 
G/N = IR has a natural C-invariant measure, which is the Lebesgue measure dz. The 
representation space of is then the completion L 2 (R, C) of the space of compactly 
supported continuous sections in the homogeneous Hermitian line bundle Ex C, and 
the action of G on it is given by 

T~ k (g)f\z] = e'^- lz ^f[{g~ l z) H } = e ^ F ^f[z-g z }, g = (g,g z ) G G, f G C (R,C), 

where for any g G G we write g = gNgH, 9n G N, gu G H . For / G Cq(G) define the 
(harmonic analytical) Fourier transform by 

/(tt) = tt(/)£) w = / f{g)7t(g)D,dg, 
Jg 

where the operator D,,- is defined on G Co (JR., C) by 

D v <f>[z] = A(z) + ^(f)[z} = (det F(z))^(j)[z]. 

(Note that there is a misprint in the formula (7.49) of [9], and the sign "— " in the 
power of A should be replaced by "+". The author of [9] confirmed this in a private 
communication.) By Theorem 7.50 in [9], the operator fields f(ir) are measurable 
fields of Hilbert-Schmidt operators, and if we identify the space of Hilbert-Schmidt 
operators on with the tensor product space H n Cg> H*, then the Fourier transform 
gives an isomorphism 

L 2 (G)~ / dv(n)H w ®H* w . 
Jg 

To find the Plancherel measure du(jr) we calculate the Fourier transforms f(Tt) directly. 
For <p G C {R,C) we have 

f{T^[r] = ! f{g)T~ k {g)D^[r]dg = 
Jg 



Harmonic analysis in Bianchi I- VII 

by a substitution g' z = r — g z 

= f, /" dg x dg y dg'J(g x ,g y ,r- g' z ) e^ F ^ (det F(r - g' z ))-\det F {g' z )Y 
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= f) / dg^M^gytr ~ {r)k ' S \det F(r - g^ 1 (det F(g' z ))- 2 <f>[g' z }. 

Jm 3 

Thus /(TV) is an integral operator with a smooth kernel 
^(r,^) = ffe[/(.,-,r-^ 

where 

)P)= / dxdyi/>(x,y)j#&. 
Jm. 2 

The Hilbert-Schmidt norm |||/(TV,)||| is given by 

\\\m)\\\ 2 = t) 2 [ drdg' z \)Cl(r,g' z )\ 2 . 
Jm. 2 

Coming back to the original variable g z = r — g' we have 



= ^ / drdg z ^ 2 [f(.,.,g z )](F X (r)k) (det F(^))" 2 (det F(r - g z )) 



$Af{,,9z)]{F x {r)k) (detF(^))- 1 (detF(r))- 1 



= h 2 / drdg z 
Jm 2 

(by Fubini's theorem) 



fl 2 / dg z {detF{g z ))- 1 / &*[/(.,., g z )\{F x {r)k) (detF(r))- 1 . 



The Plancherel formula. Now we refer to the previous section about the co- 
adjoint orbits, and note that in all cases p(k, r) = z>(fc)(det F(r)) -1 with some continuous 
non-negative function v{k) on K. We will shortly see that 

dvik) = \)- x v(k)dk (2) 

is exactly the Plancherel measure desired. Indeed, 



' dkr X m\\\f{TM k ))\\? = / dku(k) / ^(detF^))-^ 
k J k Jm. 



X dr 



3w[f(.,.,g z )}(F L (r)k>(k)) (detF(r))- 1 = 
by another application of Fubini's theorem (see [19], chapter XIII), 

= f) / ^(detF(^))- 1 [ dk f drp(k,r) $*,[/(., , ^)](F- L (r)/fe (A;)) 



K 
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by definition of p(k,r), 

= f) / ^(detF(^))- 1 / dk 

JR JR 2 

by the Euclidean Parseval formula, 

= f> / d^(detF(^)) _1 / dg x dg y \f(g x ,g y ,g z )\ 2 = / dg\f(g)\ 2 , 
Jr Jr 2 Jg 

thus we arrive at an explicit Plancherel forumla, 

fdu(k)\\\f(T k{k) )\\\ 2 = I dg\f(g)\ 2 . 
The Plancherel measures for groups II- VII are thus given by 

vn(k) = \k\, vin(k) = 1, i>iv(k) — l + k±, i>v{k) = 1, i>vi-{k) = cos 2 (k) — qsm 2 (k) 

v VI+ ( q ) = q k2 ™ d2 ,v VII (k) = \k\. 

Note that we could have chosen the cross section for VI, q < in the same way as for 
VI, q > to get a uniform Plancherel measure v V i = i / VI+ for all Bianchi VI groups, 
but we preferred the more conventional circle to the quartet of rays in Figure 1 when it 
was possible. This can be altered for any technical purposes when needed. 



5* [/(.,.,&)](*) 



6. Scalar spectral analysis on Bianchi I-VII groups 

Here the term scalar spectral analysis is understood as the spectral theory of the scalar 
Laplacian. Of course, there is no distinguished Laplacian on these groups. We will 
consider any Laplacian which arises as the metric operator with respect to any conserved 
metric on the group. 

Let G be one of these groups, and let £(G) be its Lie algebra generated by three 
right invariant vector fields £1,^2? £3- Let further X-y, X 2 , X 3 be a basis of left invariant 
vector fields on G, and du 1 , du 2 , du 3 the dual basis. Any left invariant metric h a b on G 
can be written as 

3 

hab = ^2 h ij duJ a duJ b> 

where hij is any symmetric positive definite 3x3 matrix, and the corresponding metric 
Laplacian will be 

3 

A, = h ij XiXj, (3) 

with h lj = (hij)^ 1 . To see this first note that 
3 3 3 

h ij XiXjf = ^ E \ X \ X T d ^ + x KdiX?)d m ]f. 

i,j=l i,j=l l,m=l 
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On the other hand the connection Laplacian related to the Levi-Civita connection is 
given by 

3 3 3 

A ^ = E w E \ x l x T d ^ -J2 x i x T T imdk\, 

i,j=l l,m=l k=l 

where the rf m are the Christoffel symbols. This together with the observation 

E x l x T T "m = -\i x \di x T + xl A x n 

l,m=l 

which follows from V^Xj = |[Aj,X,-], gives ([3]). Our aim will be to find the 
eigenfunctions and the spectrum of A^. If £i and £2 commute and also commute 
with all Xi, then £1,^2, A triple of commuting operators, and have common 

eigenfunctions. We will find those eigenfunctions and show that they are complete in 
the sense we desire. For the ease of notation let us denote 

h 2x2 = h lj \ij <3 , 

h' 3 = h^\ i<J=3 , Ir- h ij j , : , (4) 

First let us describe the spectrum Spec(A/ l ) of the Laplacian A^. We note that 
A h is a negative semidefinite operator, as (A h f,f) L 2^ = —(dhf,dhf)L 2 {G) < 0, where 
dh is the exterior derivative with respect to the metric h^. Thus A^ is a semibounded 
and real symmetric operator on L 2 (G). There are several ways of extending A^ to a 
self-adjoint operator on L 2 (G). A real symmetric operator has a self-adjoint extension 
by von Neumann's theorem [20] ■ A semibounded symmetric operator has a self-adjoint 
by Friedrich's extension theorem [20]. But we have something stronger. The Lie group 
G with its left invariant Riemannian metric h a b is a complete Riemannian manifold 
\21\ . Then following [22] A^ is essentially self-adjoint on C^°(G). Being a negative 
self-adjoint operator A h has a real non-positive spectrum, Spec(A /l ) C (— 00; 0]. The 
semidirect structure of our groups satisfies the conditions of Lemma 5.6 of [21], and we 
have for the scalar curvature Rh the following formula, 

R h = -Tr[S 2 } - (Tr[S]) 2 , 

where we took into account that the normal Lie subgroup M. 2 with the induced metric 
is flat. The matrix S is given by 

1 1 

S = -(ad (0A1) | R 2 + od( 0j0l i)|Ja) = 9 W) + /W*)' 

where the adjoint * is understood as 

h(Af,g) = h(f,A*g), Vf,g G £(C), A e Aut(£(G)). 

Here h(f, g) for f,g£ £{G) means the evaluation of the Riemannian metric h on the 
vector fields f,g. Thus all our groups endowed with any left invariant Riemannian metric 
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are spaces of constant negative curvature equal to Rh, which is given explicitly in terms of 
the matrices f(l) = M and h 2x2 with S = |(M + (h 2x2 )~ 1 Mh 2x2 ). This in turn implies, 
following [23], that the essential spectrum of is precisely EssSpec(A/ l ) = (—00; Rh]. 
Recall that the essential spectrum of a self adjoint operator consists of eigenvalues of 
infinite multiplicity (see [20] )■ For the group Bianchi I, all irreps are 1-dimensional, and 
as we will see later in the section, each eigenspace representation includes an infinite 
number of them, thus there is no discrete spectrum. For the remaining groups, we 
have seen in the previous section that no finite dimensional representation enters the 
Plancherel formula. On the other hand, in the next section we will see that the infinite 
dimensional eigenspaces exhaust L 2 (G), hence no finite dimensional eigenspace exists, 
i.e., the discrete spectrum is empty, and therefore Spec(A ft ) = EssSpec(A ft ). 

To find the generators & for Bianchi I- VII groups we differentiate the left translation 
map x 1 — y gx, 

g(x,y,z) = ({g x ,g y ) + F(g z )(x,y), g z + z), 

and obtain 







( 1 °\ 




(d x \ 


6 




1 






W 




[(x, y)F T (0) I) 




w 



We see that £1 = d x and £2 = d y do indeed commute. To find the left invariant vectors 
Xi (which are the generators of right translations) we differentiate the right translation 
map x 1 — y xg, 

(x,y,z)g = ((x,y) + F(z)(g x ,g y ),z + g K ), 

and get 



(xA 


(F T {z) 0\ 






x 2 









\X 3 ) 


V 1) 




w 



Thus ^1, ^2 do commute with all Xi. Now let ((x) G C°°(G) be a joint eigenfunction for 
{£i,£ 2 , A/j}. Then it is necessarily of the form 

((f) = e ifo} P(z), 

where kc E C 2 , x = (x, y), and satisfies 

A h ((x) = AC(x), 

for some A G C. A matrix representation of equation ([3]) and a bit of manipulation 
yields the following equation 

h 33 P(z) +t(klF(z)[h' 3 + (h 3 'y])P(z)- 

-(A + klF{z)h 2x2 F T {z)k c - ih 3 'F T {z)M T k c )P{z) = 0, 
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where F T (z) = d z e zM = F T (z)M T was used. This is a generalized time- dependent 
harmonic oscillator equation, which always have solutions, and those solutions comprise 
a two complex dimensional space. For given A and kc let us choose two linearly 
independent solutions P\^ c (z) and Qx,kc( z ) (the choice of initial data may be arbitrary). 

First we consider the group Bianchi I. Here M — 0, F(z) = 1 and the equation 
becomes 

h 33 P(z) + i(k^[h' 3 + (h 3 ') T ])P(z) - (A + k^h 2x2 k c )P(z) = 0. 
One can easily check that P(z) = e lkz ' z is a solution if 

A = -klh lj k c , 

where kc = (kc,k z ). This is a consequence of the fact that for this group £3 also 
commutes with all ^ and JQ, so that there exist joint eigenfunctions of the commuting 
operators £1, £2, £3, &-h of the form 

corresponding to the eigenvalues 

A = -klh ij k c . 

In particular, when we restrict ourselves to the irreducibles kc — k G M. , we obtain 

and we observe immediately that each eigenspace corresponding to the eigenvalue A 
includes infinitely many k which satisfy 

A = -k T h lj k. 



k 



Of course, e z ' z do not exhaust all solutions P(z). But it turns out that the 
constructed in this way are already complete in L 2 (G). Indeed, that is the essence 
of the Euclidean Parseval equality. To be more precise, we need to take du(k) = ^dk as 
the Plancherel measure for the Euclidean Plancherel formula to hold. Equivalently we 
can renormalize by taking 

so that the Plancherel measure is independent of the metric. But this ease of 
construction is a peculiarity which the remaining groups Bianchi II- VII do not share, 
and we proceed to determine their eigenfunctions. 

For the groups II- VII let us now restrict to > A G K and k c = F (r)k (—k) e 
M 2 \ V°, k E r E~M, (minus sign for convenience). The equation now becomes 

h 33 P(z) + i(k (-k) T F(z - r)[h' 3 + (h 3 ') T })P(z) - 
(A + k (-k) T F(z - r)h 2x2 F T (z - r)k (-k) 

-ih 3 'F T (z - r)M T k (-k))P(z) = 0, (6) 
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and the two independent solutions will be denoted by P\ t k, r ( z ) an d Q\,k,r(z)- If we se t 
P\,k,o(z) — P\,k{z), Q\ t k,o{z) — Q\,k(z), then a variable substitution z — r \-t z shows 
that we can choose Px,k,r(z) — P\,k{z — r), Q\ t k, r {z) — Q\,k(z — r )- Another point that 
can be noticed in equation <^ by taking the complex conjugate is that we can choose 
P\,-k{z) = P\k(z), Q\,-k(z) — Q\,k( z )- Finally we construct the eigenfunctions 

Cfc,A,r,-(5) = (det F{-r)y^ r ^- k ^P XAs (z - r), (7) 

where to s — 1 (— 1) corresponds P\,k,s = P\,k (Q\,k)- Note that each (k,\, r ,s enters 
with its conjugate, Cfc,V,s — C-fe,A,r,s- As we will see in the next section, P\,k,s are 
orthogonal with respect to the weight detF(— z), which shows that (k,\,r,s j us ^ defined 
are orthogonal with respect to the same weight. Again, instead of using the Plancherel 
measure (T5]) we can use dv(k) = v{k)dk and renormalize according to 

Ck,x,rM = 4F(detF(-r))e^ ± (^^(" fc )' £ >P A , fciS (^ - r). 



Note that by (jSJ) the number () is just J det hij. 



7. Fourier transform on Bianchi II-VII groups 



As a first step on the way of establishing the completeness of {(k,\,r,s} we prove a simple 
proposition. Consider the differential operator 

£>- = h 33 ^ + i(k T F(z){h' 3 + {h 3 ') T })^- - (k T F{z)h 2x2 F T (z)k - ih 3 'F T (z)M T k), 

k e M 2 \ V°, 

which by definition satisfies 



Proposition 1 The operator with domain Cq°(M.) is symmetric in L 2 (M, det F(-z)dz), 
for any k G M 2 \ V° . 

Proof: Let us first write Green's identity for the operator on the infinite tube 
D 1 x R c G where D l is the unit disk in the x-plane, 



dx (e- i{ ^ } g{z)A h \e i{k * } f(z 



D 1 xl 



e-^g(z) 



J,{k,x} 



d 



dzdl(x)(e-^g(z)(x,^) 



i{k,x} 



d 



(X. 



' dx 

dzdl(x)2ig(z)f(z)(x, k) = 0. 



S 1 xl 
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Next we note that 

Jd i xr v l j l j / 

dxdydz{detF(-z)) (g{z)D % f(z) - D k [g{z))f{z)) = 

D 1 xR 

= vr / dz(detF(-z)) (g(z)D h f(z)-D k [g(z)]f(z)) = 0, 
Jr 

which holds on the dense subset of all f,g in C£°(R) inside L 2 (R, det F(— z)dz), and 
symmetry is thus proven. □ 

Now from the definition it is clear that D k is a negative definite operator (because 
Ah is such), and is hence upper semibounded, and has a self-adjoint extension in 
L 2 (R, det F(— z)dz) by Friedrichs extension theorem [2D]- In particular, for k = ko(— k), 
k G the generalized eigenfunctions {P\,k,s} \eSp(A h ),s=±i are complete and give rise to 
a Fourier transform S^/m on L 2 (IR, det F(-z)dz) by means of an abstract eigenfunction 
expansion. 3" ko (_ k) is given by 

(ffjS (- fc) /)(A,a)= / efe(detF(-*))P^)/(4 
Define now the linear isomorphism 23 : L 2 (M) — > L 2 (R, det F(— z)dz) by 

/(z) = [QJ0](z) = 0(-z)(detF(z))i 
This induces a Fourier transform Jfc = -5'k (-ky® which acts as 

(3W>)(A,s)= / ^(detF(z))5P Aifc> .(-z)0(z) = ^(A,A;,s). 
The inversion formula is given by 

00) = (detF(»)i V / d\4>(\,k,s)P x 



\,kA~ z )- 



Now we are in the position to show how Cfc,A,r,s are related to the irreducible 
representations T^ r k y Consider the following transformation on / 6 C£°(G), 

f(k,\,r,s)= / dgCk,x,r,s{9)f(g), 

JG 

with eigenfunctions Cfc,A,r,s defined in Section [61 We will see that f(k, A,r, s) are in 
some sense proportional to the matrix columns of the operators f{T ko ^). First we see 
that 

f(k, A, r, s) = f)(det F(-r)) / dxdydz(det F(-z)) x 
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xf(x,y,z)e« F± ^ k ^P Xtk , s (-(r - z)) = 
= f)(detF(-r)) / dxdydz(det F(-z))(det F(r - z))~*f(x, y, z)e 
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i(F ± (r)k (k),x) 



X 



x(det F(r - z))*P X)k;S (-(r - z)). 

Next we recognize that this is related to the extension of the operator /(Tj^^O from 
L 2 (R) to C c 



f(k,X,r,s) = (detF(-r))/(r fe(fc) ) [(det F(r))^P AiM (-r) 
Integrating we obtain 



^ / dXf(k,X,r,s)4>(X,k,s) = (det F(-r))/(T^ o(fc) )0[r]. 

s=±1 ./Sp(A h ) 



(8) 



Recall now the Fourier inversion formula as given in [12] (notation there is different, 
and we have adapted them to ours adopted from [9]), 

(9) 



f(l) = / dv{k)Tr DJ(T h{k) ) 
J a 1 

Formally a matrix element of D v f{T kQ ^) would be an expression 

(detF(2))ip ViM ,(-«),D w /(T fe(fc) )(detF( 2 ))ip AiM (-^ 



L 2 ( 



dz(det F(z))P X ',k,s>(-z)f(k, A, z, s), 
which does not make sense in precise terms. However, the trace of such elements, 

V / dX I dz{detF{z))P x , k>s {-z)f{k,X,z,s), 

s=±1 JSp(A h ) JR 

can be given an exact sense if we change the order of integration, 

dzJ2 [ dX(detF(z))P x , k>s (-z)f(k,X,z,s). 

s =±l JSp(A h ) 

Indeed, let {p n (z)} be an orthonormal system in L 2 (R). Consider the Fourier 
transforms p n (X, k, s), and consider the following bi-distribution in the Fourier space, 
Y^=iPn{X,k,s)p n (X',k,s'). Let f,g be the Fourier transforms of arbitrary /, g G L 2 ' 
We have 

EEE / dX f dX%(X, k, s)p n (X', k, s')f(X, k, sfg(X', k, s') = 

n=1 s=±1 s , =±1 JSp(A h ) JSp(A h ) 



E E dXp n (X,k,s)f(X,k,s)) V / dX'p n (X',k,sTg(X\k,s') 
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oo 
n=l 

thus Y,n=iJn{\ k, s)p n (X', k, s') = 6(X - X')5 S S ,. Now 
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/ dXf(X,k,s)g(X,k,s), 

5=±1 JSp(A h ) 



dzY^ / d\(detF(z))P x , k!S (-z)f(k,\,z,s) 

s=±1 JSp(A h ) 



oo p p 

dz(detF(z))J2J2 / / d\d\%(\,k,s)p n (\',k,s')P XAs (- 

n=l s,s' J , . so 



z)f(k,X',z,s') 





OO 


( 






I dz{detF{z))J2 


^2 / dXp n (X, k, s)P Xtk)S (-z) 


X 




n=l 


\ S=±1 S P (A h ) 


J 






( 


\ 




X 


E / 


d\'f(k, A', z, s')p n (X', k, s') 






\ s ' =±1 Sp<[A h ) 


J 





using USD, 

r oo 

dz(det F{z)) £(det F(z))-^(z)(det F(-z))/(T^ o(fe) )p n (z) = 

n=l 

as both the sum and the integral converge in L 2 , 

oo „ oo „ 

n=l ^ 

z)f(k,X,z,s). 



n=l 



= ^{Pn,D^f{F kQ{k) )p n ) L 2 m = Tr 

71=1 

Hence from ([9]) we have 

f(l) = [ dv(k) [ dzJ2 [ dX(detF(z))P XAs (- 

J& JR s=±1 JSp(A h ) 



To find an inversion formula at an arbitrary point g G G we apply this to the left 
translated function [L g -if]([x) = f([gx), 

f{g) = \L g -,f){l) = [ dv{k) [ dz V / dX(detF(z))P XAs {-z)[L^Tf){X,k,z,s). 
In J® s=± i Js P (A h ) 

But from the definition 



[L g - 1 f](X,k,r,s)= / ^Cm^(^)[V 1 /](^)= / dh , Ck,x, r ,.(9- 1 ti)f(h'). 
Ig Jg 
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From the definition of (k,\,r,s we find 

Cm^GtW) = e -^>+^°«'*)(det F(g z ))(xAr +9 Uh'), 

thus 

' dh'Ckx^g-^fih') = e^ F± t r+ °'foW*\detF(g z ))f(k,\,r + g z ,s). 



G 



Therefore 



f(g) = [ dv{k) f dzJ2 I dX{detF(z))P Xjk>s {- 

J& JM s=±1 JSp(A h ) 



z)x 



xe ^(F^(z+g z )k (k),g)( det F(g z ))f(k, X,Z + 9z , s) 

by substitution r = z + g z 



dv{k) [ drY] I d\(detF( 
s=±1 JS P (A h ) 



r))f(k,X,r,s)x 



dv{k) [ drJ2 [ d\{detF{r))f{k,\,r,s)( kAr>s {g), 

Jm s=±1 JSp(A h ) 

which is our final inversion formula. 

It remains to note that by denoting a = (k, X,r, s) we have satisfied all conditions 
for the eigenfunction expansion Ca{f) to give a conventional Fourier transform in sense 
of 0. 

8. Automorphism groups of Bianchi I-VII groups 

In this section we consider the automorphism groups Aut(G) of Bianchi I-VII groups. 
After performing the calculations we discovered that these automorphisms have been 
obtained earlier in [21]. However we give here also the dual actions of these 
automorphisms on G which is new. This may become important when analyzing the 
transformation in the Fourier space induced by automorphisms. We start by noting that 
Bianchi I-VII groups are matrix groups, and their matrix realization can be given by 



y 

W 



H> G(x,y,z) 



( x\ 

F(z) 

y 

o o 1 



It can be easily seen that in this realization the group multiplication indeed corresponds 
to the matrix multiplication. The respective Lie algebra realization will be 



y 

W 



H> Q(x,y,z) 



zM 



y 



\o o oy 
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which again can be checked to intertwine the matrix commutation with the Lie 
bracket. Moreover, we could have obtained immediately the exponential map by setting 
exp(x,y,z) = exp(g(x, y, zj) instead of referring to the Zassenhaus formula, but the 
latter is a more Lie theoretical approach. Now that all Bianchi groups are connected 
and simply connected by Theorem 1 of III. 6.1 in [25] it follows that Aut(G) = Aut(g) in 
the sense of a topological group isomorphism (see also [15]). An algebra homomorphism 
of matrix algebras is necessarily linear in the matrix elements. It follows that any 
at G Aut($j) depends linearly on x, y, z, and is therefore given by an affine transformation 
in R 3 , which is actually a linear transformation because it preserves 0. Therefore we 
first determine Aut(g). Let the linear map a : R 3 — > R 3 be given by 











_ ( «2x2 a, 3 \ 


y 


\ "3. "33 J 


r 


V) 




W 



where we use notation similar to (J4j). Then a G Aut(g) if and only if a[x, y\ = [ax, ay\, 
where [, ] is the Lie bracket. Expanding this condition we get the system of requirements 

"2x2^ - a 33 Ma 2x2 + Ma, 3 a 3 . = 0, (10) 

a 2x2 MaaJ, = 0, 

a 3 .M = 0, 

where a is the unit antisymmetric matrix. The patterns of admissible matrices a 
satisfying this system have to be computed for each group independently. For Bianchi 
I we have M = and all three conditions are satisfied trivially. For Bianchi IV- VII the 
matrix M is invertible hence the third requirement means a 3 , = 0, so that the second 
becomes trivial, and the first reduces to a 2X 2M — a 33 Ma 2x2 = 0. The cases of groups 
Bianchi II and III are a bit more involved, but the calculations are straightforward. We 
present the results in the Table El Note that whenever a 3 , = the invertibility of a 
requires a 33 ^ 0. As it can be seen from the table some algebras allow for reflective 
automorphisms and their automorphism groups consist of two components (this is what 
the union symbol |J in Table [8] refers to). Matrices of these pattern forms exhaust the 
groups Aut(g). One can compare this pattern of automorphisms to those available in 
the literature, for instance, of the Heisenberg algebra in [26J. 

Now the corresponding group homomorphisms A G Aut(G) can be found by 
composing a G Aut(g) with the exponential map, Aexp((x, y, z)) = exp(a(x, y, z)). 
Recall that the exponential map is given by 

exp((x,y,z)) = ([1 + F(z)D(z)](x,y),z), 

(where D(z) is defined in ((T]))and because this map is bijective we know that the matrix 
[1 + F(z)D(z)] is invertible for all z. The logarithmic map can be written as 

log((x, y, z)) = ([1 + F(z)D(z)]~ 1 (x, y),z), 
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Table 3. Patterns of permissible matrices a for Bianchi I- VII algebras 



and the action of the group homomorphism A related to the algebra homomorphism a 
becomes 



.4 



x 



[1 + F{z')D{z')\ (a 2x2 [l + F(z)D(z)}- 1 x + a. 3 z) 



z' 



z' = a 3 .[l + F{z)D{z)\- l x + a 33 z, 
for Bianchi II- VII groups and 



A\ = a 

z \ z 



for the Bianchi I group. From 6t 3 ,M = it follows that a 3 ,[l + F(z)D(z)] 1 = 03.. 
Thus the formula for Bianchi II- VII simplifies to 



A 



x\ _ f[l + F(a 3 ,x + a 33 z)D(a 3 .x + a 33 z)} (a 2x2 [l + F{z)D{z)\ 1 x + a, 3 z) 

a 3 .x + a 33 z 



One more step can be done in this generality. From formula ( flO|) and a 3 , M = it 
follows that 

a 2x2 M m = (a 33 M) m a 2x2 
for m > 2 and therefore for any sequence of complex numbers {fm\m=o 

00 00 
d 2 x 2 f™ Mm = J2 /-(«3 3 M) m d 2x2 + hMa. 3 a 3 . 

m=0 m=0 

whenever the left hand side exists. This can be used to establish that 

[1 + F(a 3 ,x + a 33 z)D(a 3 ,x + a 33 z)}a 2x2 = a 2x2 [l + F(^— + z)D(-^— + z)]. 

«33 "33 

This far on the explicit form of the group automorphisms. 

Now let us look at the dual spaces G. If A £ Aut(G) and 7r e G then tt o A = it' for 
some 7r' G G. Thus A induces a pullback map A* : G — > G. Because dim7r = dim7r' it 
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follows that A* maps generic representations into generic representations and singletons 
into singletons. The representations it G G are in a bijective correspondence with the 
derived representations dir which are irreducible representations of the Lie algebra q. 
In a similar fashion, any a G Aut(g) induces a pullback map a* : dG — > dG between 
derived representations. This pullback map is easier to study than that for the group 
representations. Consider first the Bianchi I group. The irreducibles are given by 

T % {g) = 

and the derived representations are 

dT^(x) = i{k,x}. 

An automorphism x = aq induces the pullback map a*(k) = a T k. Consider now the 
singletons of a Bianchi II- VII group. They are given for k G V° © M by 

T^g) = e i{ ^ } = e i{fc ' s} e ik39 % 

and the derived singletons are 

dT^(x) = i{k,x}, 

and again, an automorphism x = a q induces the pullback map a*(k) = a T k. This in 
particular means that k! = cij x2 k + ^aj,, and if k G V° then 

M T k' = M T a] x2 k + k 3 M T aJ. = 0, 

where (|T0|) and a 3 ,M = were used. We explicitly observe that the automorphisms map 
singletons into singletons, as expected. Finally we turn to the generic representations. 
Let be a generic representation of G. Then it acts on L 2 (IZ) by 

T~ k (g)f[w] = e^ F ^f[w - g z ], g = (g,g z ) G G. 

Its derived representation will be 

dT~ k (x)f{w] = i{k,F(-w)x}f[w] - zd w f[w}. 

Under the automorphism x = aq it will turn into 

dT k (q)f[w] = i{k, F(-w)[a 2X 2q + a. 3 s]}f[w] - [a 3 .q + a 33 s}d w f[w]. 

For simplicity we will consider only the automorphisms with a 3m = 0. Thus we omit 
only some automorphisms of the Heisenberg group, but this group is a central subject 
in harmonic analysis, and the missing results can be found in the literature. Define the 
isometric isomorphism X : L 2 (R) — > L 2 (M.) by 



'ot 33 
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Consider the representation dTy with k' = c\ 2 T x2 k. Note that because = we have 
that aj x2 is invertible, and from ([TO]) we assure that it maps k ^ V° to k! ^ V°. Thus 
dTy is generic. Its action on the image %{f)[w] is given by 

dT^q)%U)H = i{k',F(-w)q}Z(f)[w] - sd w Z(f)[w] = 

= %( i {k',F(-— )q}f)[w]+%(is{k,F(-*)a. 3 }f)[w}-Z(sa 33 df). 

«33 

Recall that we have seen that from ( ITUl) and a 3 , = it follows a 2 x2F(z) = F(a 33 z)a 2x2 , 
hence 

i{k', F(-—)q} = i{a T 2x2 k, F(-—)q} = i{k, F(-.)a 2x2 q}. 

"33 «33 

We finally see that 

dT~ k ,(q)Z(f)[w] = % (\i{k,F(-.)a 2x2 q} + is{k, F(— )&. 3 }}f - sa 33 df^) =%{alT^q)f) , 

which means that T intertwines the irreducible representations dT^oa and dT &T ^. Thus 

these two representations are unitarily equivalent, a*(k) = aJx 2 ^- If the cross sections 
are chosen explicitly (for instance, as we did) then it is a straightforward calculation 
to find the action of a* on K and ^. We omit these calculations here because, first, 
they depend on the preferred choice of the cross sections, and second, they involve 
transcendental functions (e.g., the solution of the equation e y + ay = x) and are not 
transparent visually, and do not provide a better insight into the matter. 

9. Separation of time variable in homogeneous universes 

We want to see to which extent the technique of mode decomposition developed in [6] is 
applicable to hyperbolic fields on Bianchi type and FRW cosmological models. For this 
aim we have to check whether the conditions of Proposition 2.3 are satisfied. Recall 
that the metric g of a homogeneous spacetime M = X x S, where X is an open interval 
and E is a Bianchi type homogeneous space, is given by 

3 

ds 2 g = dt 2 — ^ h a/ 3(t)du a (x)duj l3 (x), 

a,P=l 

where hij{t) (t e X) is a smooth positive definite symmetric matrix function, and did 1 
are the left invariant 1-forms on S. Condition (i) of Proposition 2.3 is automatically 
satisfied because goo = 1- For condition (ii) note that 

E n-'u^ir) E E ^(t)^(t)^(^)^(^)(^) l (f)(^) J (f) = 

i,j=l i,j=l a,/3=l 7,5=1 

3 3 

= E E h a P(t)ti 7S (t)(X a ,duj'i) h (Xp,duj s ) h = 

a,P=l 7,5=1 
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3 3 

= EE h^{t)n lS {t)5% = Tr[h-\t)kt)\. 

a,P=lj,S=l 

This shows that the condition (ii) is also satisfied. We see that the homogeneous 
spacetimes are an ideal playground for mode decomposition. Note that FRW spacetimes 
correspond to the choice h a p(t) = a 2 (t)<5 Q( g. 

If conditions (iii) and (iv) are also satisfied depends on the chosen connection V. 
For the scalar field (iii) is automatically satisfied with T = 0. In |6] we defined the field 
operator D = □ + m 2 (x), and the instantaneous field operator Dj] t = —A + m 2 (x). The 
eigenfunction equation of the latter operator was written as -Ds t C« = K(&)Ca- Condition 
(iv) of Proposition 2.3 in [6J can be obviously satisfied if h a p(t) = a 2 (t)h^ as in this 
case the time evolution amounts only to a rescaling of A (a) in D^Xa = At («)£«. Note 
that because D-£, t is G-invariant, the term m 2 is a function of t only. This is the situation 
where the dynamics of the universe consists of merely an isotropic rescaling. Thus, for 
instance, in case of FRW spacetimes condition (iv) is satisfied automatically. 

But the condition (iv) can be also satisfied non-trivially with an anistropic rescaling 
and even some shears and rotations. This is clearly possible for Bianchi I group, because 
the eigenf unctions do not depend on the matrix h. For Bianchi II- VII groups one has 
to look at the equation (EJ) to see to which extent the solution P(z) depends on the 
matrix h. Suppose h and j are two matrices for which there exist two common linearly 
independent solutions P(z) and Q(z). Because we have already seen that an isotropic 
rescaling is always possible, without loss of generality we assume h 33 = j 33 (we again 
use the notations H]). Fix k G R 2 \ V° and < A G R. Now the condition that the two 
equations have the same solution spaces can be cast into the following pair of equations, 

h 3, F T {z)k = ] 3, F T {z)q, 

A + k T F(z)h 2x2 F T (z)k = A' + q T F(z)j 2x2 F T (z)q 

for some gGM 2 \l /0 ,0<AGlR and for all z G R. That non-trivial possibilities exist is 
clear visually, but we will not go into details here. Once these conditions are satisfied 
at all t G I for the 1-parameter family of matrices h(t) describing the evolution of the 
spatial metric then condition (iv) is satisfied, and we have an explicit formula for the 
time dependent eigenvalue At (a). 

As electromagnetism is of primary importance for us, let us finally show that 
the assumption h a p(t) = a 2 (t)/i° /3 is sufficient to satisfy condition (iii) for the 1-form 
field. Indeed, the 1-form field is given by the Levi-Civita connection, for which the 
connection forms are = — T? b . Let us compute the symbol Tq 6 . It is easy to see 
that T° ob = Tg = 0. For a, b > we have 

m=l o,/3,7=l 
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If h a/3 (t) = a 2 (t)h a/3 then d t h a/3 (t) = 2H(t)h af) {t) with H(t) 
constant, and we get 

3 

r o a 6 = 2H(t) £ h^(t)h^(t)X a a (x)d^ b (x) = 2H(t)5 a b . 

a,/3,7=l 

Thus To = — 2H(t)0®l3 is not only a function of t, but also commutes with any matrix, 
hence (iii) is trivially satisfied. 

10. The mode decomposition of the Klein-Gordon field 

The investigation of classical and quantum fields in anisotropic cosmological models has 
been carried out by different authors since decades (see, e.g., [27].[28].[29].|30]). However 
these works mainly concentrate on Bianchi I models where the harmonic analysis and 
mode decomposition are obvious. With the background developed above we can extend 
this to all Bianch I- VII models. 

In [B] we have described how the mode decomposition can be performed explicitly 
for an arbitrary vector valued field given the explicit spectral theory of the model spatial 
sections E. As we have already explicitly constructed the spectral theory of the line 
bundle over Bianchi I- VII spacetimes, we can apply the mode decomposition to the 
Klein- Gordon field and see what can be gained by this technique. 

Let M be Bianchi I- VII type spacetime, i.e., a 4-dimensional smooth globally 
hyperbolic Lorentzian manifold with a smooth global time function chosen [31], and 
with the isometry group G which is one of the groups Bianchi I- VII, so that G acts 
simply transitively on the equal time hypersurfaces S t (for all missing details see [6]). 
Or to put it in simpler words, let topologically M = MxG where G is a Bianchi I- VII 
group, and let the metric be given by ds 2 = dt 2 — h^{t,x)dx l dx\ x = (x 1 , a; 2 , x 3 ) G S t , 
so that for any t £ R the Riemannian metric hij(t, .) is left invariant under the action of 
the underlying Lie group G. In [6] we cosnidered the field operator D = □ + m 2 , where 
□ is the d'Alembert operator related to the Levi-Civita connection, and m 2 e M + is a 
positive constant. The Klein-Gordon field is described by the equation 

D(p = (n + m 2 )0 = 0. 

As we have seen in the previous section for each Bianchi type there are certain 
restrictions on the dynamics of the spatial metric h(t) for the mode decomposition 
to be applicable. Recall that h(t) is described by the positive definite matrix h(t). For 
simplicity let us consider only an isotropic rescaling, 

h{t) = a 2 {t)h{0). 

We further wrote in [6] D = D t + -De ( , where D t is a differential operator in variable 
t, and _Ds t is the instantaneous field operator D^ t = —A t + m 2 . Its eigenfunctions 
satisfying D^ t (a = ^a(t)(a will be the eigenfunctions of the Laplace operator, — A t ( a = 



30 

= ?jpr being the Hubble 
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(X a (t) — m 2 )C, a . Because the time dependent Fourier transform is normalized at time 
t = 0, we have A Q (0) = — A + m 2 where a = (k,X,r,s). Using h lj (t) = a~ 2 (t)h lj (0) it 
follows that At = a~ 2 (t) Ao and hence 

A Q (t) = -^ + m 2 . 

As we have demanded that the spacetime M = X x E is globally hyperbolic, the Klein- 
Gordon field operator D = D+m 2 has unique advanced(+) and retarded(-) fundamental 
solutions E ± : C °°(M) ^ C^(M) with the properties 

(i) (□ + m 2 )E ± f = f = E ± (D + m 2 )f 

(ii) snpp{E ± f} C J^supp/) 

for all / G Cq 3 (M). Here, J^(N) denotes the causal future(+) and past(-) of a 
subset N C M. We refer to [32] for full discussion and further references. Then 
E = E + — E~ is called the causal propagator of the Klein-Gordon operator □ + m 2 
on (M,g). Any <fi = Ef, f G C^°(M) is a solution of the homogeneous Klein-Gordon 
equation (□ + m 2 )0 = 0, and the restriction of to any Cauchy surface is compactly 
supported. We define Sol (M) = EC^(M). We also write E(f,h) = (f,Eh) L2{M) 
where 

(f,h) L 2 (M) = f(x)h(x)dvol g (x); 



>M 

dvolg is the volume form on M induced by the metric g. Moreover, we set /C = 
Cq°(M)/ keiE. Then (the real part of) K becomes a symplectic space with symplectic 
form 

a([f], [h]) = E(f, h), [f} = f + ker E, [h] = h + ker E. 

The map K i-> Solo(M), [f] Ef is a symplectomorphism upon endowing SoIq(M) 
with the symplectic form 



JC 

for any Cauchy surface C in M having future-pointing unit normal field n a and metric- 
induced hypersurface measure drjc- Again, we refer to [32] for a complete discussion and 
full proofs. 

Now by Proposition 2.3 in [6] any <p G SoIq(M) can be written as 

4>{t,x) = J dfi(a) [a*(a)T a (t)( a (x) + b*(a)f a (t)C a (x)] , 

where a = (k,X,r,s) and dfi(a) = du(k)d\F(r)dr. The modes T a are to this point 
arbitrary /^-measurable solutions of the mode equation 



f a (t) + F{t)t a (t) + G a {t)T a {t) = 
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such that T a and T a are linearly independent solutions. As found in [5] for the scalar 
field on an isotropically expanding universe, 

F(t) = Pit) = 3H(t) = 3^j|, G a (t) = X a (t), I(t) = a 3 (t). 

The propagator E[f] is obviously a weak solution of the Klein-Gordon equation. If we 
want to mode decompose it we have to satisfy (2.19) of [6]. By Proposition 2.6 of 
[6] we can do it if M is an analytic manifold. But the spatial metric h is analytic in 
the spatial variable x because it is expressed in left invariant fields of the Lie group 
G. Thus we only need to choose a(t) a real analytic function. The spectrum of D^ t is 
strictly uniform with the prescription u(a) = —A. Now if we restrict the initial data 
T a (0) = p(X) and T a (0) = q(X) where p,q G *4(Ho) and pq — pq = i then Proposition 
2.9 is applicable. Suppose this is done, now from the Section 2.6 of [6] we find that 
(for the line bundle obviously s(a) = 1) 

E[f]{x) = i j dn{a) [{f a ( a , f) M T a {t)( a (x) - (T a ( a , f) M f a (t)( a (x)] . 

Now we proceed to the quantization. The mode decomposition of arbitrary CCR 
quantum fields is discussed in [13] which provides a generalization of the works by |33j, 
[34] , [35] . Here we summarize some results applied to the quantized Klein-Gordon field. 
The latter is given by the field algebra A generated by the unit 1 and the elements 4>(f) 
satisfying 

(i) ^af + h) = a ( p(f) + ( p(h), 

(ii) = <p\ 

(iii) [<P(f),<P(h)} = -t-E(f,h)l, 

(iv) 0((D + m 2 )f) = 0, V/, h G V(M), aeC. 

A state w of the field is a linear functional u G A' such that w(l) = 1 and 
lu(A*A) > for all A G A. The 2-point function u)2 of a state u is the bilinear form 
u 2 (f,h) = u((j)(f)(f)(h)). It follows that w 2 ((D + m 2 )/, h) = u 2 (f, (□ + m 2 )h) = 0, 
w 2 (/, /) > and u 2 (f, h) - u 2 {h, f) = -iE(f, h). Moreover, 

MfJ) = WUYWr) = ^mWIW) = w(^(/)) = <*2{h, /), 

bj 2 is hermitian. A quasifree state a; is a state which is completely determined by its 2- 
point function oj 2 (for precise definitions see, e.g., [36], [37]). Being a weak bi-solution of 
the field equation u 2 can be mode decomposed and that is done by Proposition 4.1 of 
[13] . (For earlier results using mode decomposition of 2-point functions see [33] , [38] , [39].) 
If we denote for convenience 

;» = (T a ( a ,f) M , r(a) = (f a ( a j) M , 



where 



(f,h) M = / dxf(x)h(x), 

J M 
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then a 2-point function can be written as 



w 2 (/, h) 



a"(f u , h u ) + a»(T, h) + b»(f u , h v ) + b-(f, h) + $(P, h v ), 



where a w and W are bi-distributions satisfying certain symmetry and positivity 
conditions, and 5 is the bi-distribution given by the integral kernel of the usual delta 
function. But it differs from the delta function because we identify kernels with bi- 
distribution using the pairing 



This strange convention is chosen only for calculational purposes and can be transformed 
to the usual form if needed. In particular, the 5 term in the formula for the 2-point 
function is 



This term by itself represents a quasifree pure state, and the remaining part of the 
generic 2-point function is symmetric. This state depends on the choice of the modes 
T a . Choosing different modes S a we will find a rich supply of such pure states. Then 
we can transform these states back to our original T a as follows. Let the old and new 
modes be related by S a = fi a T a + u a T a with \fi a \ 2 — \v a \ 2 = 1. The pure state given 
by the S term in modes S a is determined by the choice a u = b u = 0. As described in 
[T3] in the original modes these components will become a u (a,^) = 5(a, f3)fi a i , a and 
= 5(a, f3)\i , a \ 2 . But such states do not exhaust all pure quasifree states. By 
Corollary 4.1 of [13J any pure quasifree state is given by a w (ct,/3) = —5(a,f3) o S v,u 
and = —6(a,p) ° S v,v , where the linear maps S v ' u , S v ' v are subject to certain 

conditions. In fact, the pure states given by 5 terms are special in that they are 
homogeneous states, i.e., they are invariant under the isometry group G. More generally, 
by Proposition 4.3 of p3] any homogeneous quasifree state is given by coefficients a u , 
b u which are of the form 



with distributions a w {k, A, A', s, s') and b w (/c, A, A', s, s'). 

Another important notion is the notion of Hadamard states. A quasifree state u is 
said to be Hadamard if its 2-point function 0J2 satisfies the microlocal spectral condition 
(fiSC) (see [40J). Hadamard states are believed to be the states of physical importance 
for several reasons discussed in the literature ( [36], [H], [31], [32]). One way of checking 
whether a given 0J2 satisfies the 11SC is to try to compute its wave front set directly [40J. 
The mode decomposition suggests another way of doing this. By Proposition 4.6 of 
[T3] there exists a wide variety of modes T Q such that the homogeneous pure states given 
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by 5 terms in these modes are Hadamard. Such modes can be computed easily. For some 
conceptually related approaches (however emphasizing somewhat different aspects), see 
also [33] , [13] , [B] , [38] , [39] . For convenience we switch to the variable 

as described in [13J, then the mode equations become 

f a (s) + A a (s)T a (s) = 0, 

where 

A a (s) = a 6 (t(s))( a2 ^ +m 2 ), 

and t(s) is the inverse function of s(t). Now choose an arbitrary non-negative function 
r\ G Co°[— 1,0] and let fj(s) = f* dar}(cr). Denote 

A' a (s) = (l-fj(s))A a (-l)+fj(s)A a ( S ), 

and let 

f^s)+A' a (s)T^s) = 0. 

Choose the initial data T' a {— 1) = p(A a (— 1)) and T£(0) = q(A a (— 1)) with corresponding 
functions p, g as before with the additional constraint that 

P(A) - -== = o(A-«), g(A) - = o(A—), 

where = o(A~°°) means of rapid decay at A — > +oo. Then the choice T a (0) = T^(0) and 
T a (0) = T^(0) yields modes T a such that the homogeneous pure state given by the S 
term is Hadamard. This is basically the essence of Proposition 4.6. Some relations of 
the rapid decay in the Fourier space with the Hadamard property can be found in [38] . 
It follows that if such modes T a are chosen, then a generic quasifree state is Hadamard if 
and only if the remaining symmetric part of UJ2 given by coefficients a u and is smooth. 
A general criterion for smoothness of a distribution in terms of it Fourier transform is 
unfortunately not known in harmonic analysis. In principle Proposition 4.4 of [B] 
could give at least sufficient conditions, but this is still a work in progress. 
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Appendix A. Induced representations 

Although the generalities of induced representations can be found in any standard 
textbook on group representations, for consistency we will very briefly give an overview 
of them in this appendix. We will follow the Chapter 6 of [9] in our treatment. 

The inducing procedure produces unitary representations of a locally compact 
topological group G from a unitary representation of its closed subgroup H. In this work 
we will be interested only in Lie groups, so that the majority of functional analytical 
issues are automatically settled. If p is a unitary representation of H, and v the unitary 
representation of G induced from H (described below), we will write v = Ind^p. In 
particular the restriction of v to H is unitarily equivalent to p, u\h ~ P- Unless the 
homogeneous space G/H is a finite set, v is infinite dimensional. 

Denote M = G/H. We present the construction of the induced representation 
under assumptions which hold in cases of our interest. Namely, we suppose that there 
exists a G-invariant measure dxM on M. Any x G G can be uniquely written as 
x = xm^h with xm G M and xh G H. If H. p is the representation Hilbert space of 
p, then the representation Hilbert space H, v of the induced representation is taken to 
be H, v = L 2 (M,T-L p ,dxM), i-e., % p - valued cixAf- square integrable functions on M. The 
action of the representation v on 7-L v is given by 

v{x)f{y) = p{{x- x y) H )- 1 f{{x- 1 y) M ), VxeG,yeM, f G U v . 

That this is a natural construction can be seen by the following nice properties. If 
p and p' are unitarily equivalent unitary representations of H, then v = Ind^p and 
v' = Ind^p' are unitarily equivalent unitary representations of G. Moreover, it can 
be shown that if {pi} is a family of unitary representations of H then Ind^ (f) p { is 
unitarily equivalent to @Ind^pj. 
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